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ABSTRACT 

The origin of huge infrared luminosities of ultraluminous infrared galaxies 
(ULIGs) is still in question. Recently, Genzel et al. made mid-infrared (MIR) 
spectroscopy of a large number of ULIGs and found that the major energy 
source in them is massive stars formed in the recent starburst activity; i.e., ~ 
70% - 80% of the sample are predominantly powered by the starburst. However, 
it is known that previous optical spectroscopic observations showed that the 
majority of ULIGs are classified as Seyferts or LINERs (low-ionization nuclear 
emission-line regions). In order to reconcile this difference, we compare types of 
emission-line activity for a sample of ULIGs which have been observed in both 
optical and MIR. We confirm the results of previous studies that the majority 
of ULIGs classified as LINERs based on the optical emission-line diagnostics 
turn to be starburst-dominated galaxies based on the MIR ones. Since the MIR 
spectroscopy can probe more heavily-reddened, inner parts of the ULIGs, it 
is quite unlikely that the inner parts are powered by the starburst while the 
outer parts are powered by non-stellar ionization sources. The most probable 
resolution of this dilemma is that the optical emission-line nebulae with the 
LINER properties are powered predominantly by shock heating driven by the 
superwind activity; i.e., a blast wave driven by a collective effect of a large 
number of supernovae in the central region of galaxy mergers. 

Subject headings: galaxies: emission lines - galaxies: starburst - galaxies: active 
- galaxies: Seyfert - infrared - emission 
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1. INTRODUCTION 

The origin of huge infrared luminosities of ultraluminous infrared galaxies (ULIGs) 
is still in question; i.e., which is the major energy source, starbursts or active galactic 
nuclei (AGN) ? (Sanders et al. 1988a; Condon et al. 1991; Majewski et al. 1993; Smith et 
al. 1998; see for a review Sanders & Mirabel 1996). Since their far-infrared luminosities, 
L(FIR) ^ 1O 12 L , are comparable to those of quasars, possible starburst-AGN connections 
have also been discussed by many authors (Sanders et al. 1988a, 1988b; Norman & Scoville 
1988; Terlevich et al. 1992; Lonsdale, Smith, & Lonsdale 1993; Lonsdale et al. 1998; see, 
however, also Section 4 in Taniguchi 1997). 

Optical spectroscopy has been often used to investigate the activity in emission-line 
galaxies (e.g., Veilleux & Osterbrock 1987). In fact, optical spectroscopic observations 
have shown that the majority of the ULIGs are AGN-like objects; Seyferts or LINERsQ 
(Sanders et al. 1988a; Armus et al. 1989; Kim et al. 1995; Kim, Veilleux, & Sanders 1998; 
Veilleux et al. 1995; Veilleux 1997). Further, in order to find hidden broad emission line 
regions (BLRs), which provide definite evidence for AGN, near-infrared (NIR) spectroscopic 
programs have also been conducted (Goldader et al. 1995, 1997a, 1997b; Veilleux, Sanders, 
& Kim 1997). However, hidden BLRs have been found in several ULIGs. Therefore, the 
majority of AGN-like ULIGs are type 2 Seyferts (S2s) or LINERs although the fraction of 
type 1 Seyferts (Sis) increases with increasing L(FIR) (Kim et al. 1998). ULIGs which are 
optically classified as Sis or S2s are not necessarily powered by an AGN. Their classification 
only indicates the existence of an AGN. This AGN may not dominate the total bolometric 
luminosity. 

Furthermore, since some LINER-like ULIGs show the extranuclear LINER emission 



LINER = Low Ionization Nuclear Emission-line Region (Heckman 1980). 
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(Veilleux et al. 1995; Kim et al. 1998; Heckman et al. 1996; Taniguchi & Ohyama 1998), 
the origin of LINER-like ULIGs is controversial because this emission can be attributed to 
shocks caused by the interaction between starburst-driven outflows and the ambient gas. It 
is also observed that a number of ULIGs have intense star forming regions in their central 
regions (e.g., Shaya et al. 1994; Surace et al. 1998). For example, the main energy source in 
Arp 220, one of the archetypical ULIGs, is attributed to the intense starburst in the central 
region (Shaya et al. 1994; Larkin et al. 1995; Skinner et al. 1997; Iwasawa 1998; Scoville 
et al. 1998; Taniguchi, Trentham, & Shioya 1998; Shioya, Taniguchi, & Trentham 1998). 
However, possible evidence for the hidden AGN in Arp 220 has also been discussed because 
of the presence of compact OH megamaser sources in its central region (Diamond et al. 
1989; Lonsdale et al. 1994, 1998). This raises another important question; what powers 
the ULIGs on average ? (Genzel et al. 1998). Although this is actually an important 
question, it is difficult to give a certain answer because all the ULIGs have a large amount 
of molecular gas and dust in their central regions (Scoville et al. 1991; Scoville, Yun, & 
Bryant 1997; Downes & Solomon 1998 and references therein) and thus the putative central 
engine could be obscured even if present. Therefore, if we would like to probe the very inner 
regions of ULIGs, we have to perform MIR and FIR spectroscopy of ULIGs. 

Recently, Genzel et al. (1998) made MIR spectroscopy of a large number of ULIGs. 
Since the MIR spectroscopy can probe heavily reddened (A v ~ 50 mag) regions, it would be 
very powerful in searching for hidden AGN in ULIGs (see also Lutz et al. 1996). However, 
they found that their major energy source is massive stars formed in the recent starburst 
activity; i.e., ~ 70% - 80% of the sample are predominantly powered by the starburst. On 
the other hand, as described before, the previous optical spectroscopic observations showed 
that the majority of ULIGs are classified as Seyferts or LINERs. In order to reconcile this 
difference, we compare types of emission-line activity for a sample of ULIGs which have 
been observed in both optical and MIR and discuss a possible resolution of this dilemma in 
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this paper. 

2. A BRIEF SUMMARY OF THE MIR SPECTROSCOPY BY Genzel et al. 

(1998) 

Genzel et al. (1998) performed the MIR spectroscopic survey of 15 ULIGs using both 
short wavelength grating spectrometer (SWS) and ISOPHOT-S spectrometer on board 
Infrared Space Observatory (ISO). They utilized either [Ne v]A14.3/xm/[Ne n]A12.8/xm or 
[0 lv]A25.9/im/[Ne ll]A12.8//m emission-line intensity ratio as the discriminator between 
starbursts and AGNf]. Since the poverty of high energy photons from massive OB stars, 
high-ionization lines such as [Ne v] are generally weak in starburst galaxies. On the other 
hand, such high-ionization lines are often strong in AGN because of a larger number of 
high-ionization photons (Spinoglio & Malkan 1992; Voit 1992a). In fact, Genzel et al. (1998) 
obtained [Ne v]A14.3/xm/[Ne ll]A12.8 / um < 0.01 for a sample of nearby typical starburst 
galaxies (e.g., M82) while ~ 0.1 - 1 for a sample of nearby AGN (e.g., NGC 4151). 

They detected the [Ne v] emission only from Mrk 273 among the 15 ULIGs (see Table 
1). Its [Ne v]/[Ne n] ratio, 0.27, is consistent with the observed range for AGN. Although 
the [O iv] was detected in IRAS 23128-5919, its [O iv]/[Ne n] ratio, 0.037, is so small 
that this ULIG is a starburst-dominated galaxy. Upper limits for the remaining objects lie 
in the intermediate range between the starbursts and AGN, indicating the intermediate 
nature of the ULIGs. In order to estimate the relative contribution of starbursts to the 
total energy, they also used the 7.7 /im emission feature which is one of unidentified IR 

2 They note that the [S lll]A33.5;um line can also be used instead of [Ne ll]A12.8//m 
whenever the [Ne n] line data are not available because the [S ill] emission arises nearly 
the same ionization-condition region as that of the [Ne n] emission. 
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emission bands, probably emitted by polycyclic aromatic hydrocarbons (PAHs). Since these 
features can be seen only in star-forming galaxies, they are useful in estimating the star 
formation activity in galaxies (Moorwood & Oliva 1988; Mouri et al. 1990; Voit 1992b). 
Then using the diagram of the [O lv]/[Ne n] intensity ratio vs. the relative strength of 7.7 
(Km PAH feature, they concluded that 70% - 80% of their ULIG sample are dominated by 
the starburst component. 

3. OPTICAL vs. MIR CLASSIFICATIONS OF ULIGs 

We compare the optical emission-line classification with the MIR one for a sample of 
ULIGs which are observed in both optical and MIR. Our sample consists of 12 ULIGs given 
in Table 1. We have compiled the following important MIR and optical emission-line ratios; 
a) MIR: 1) [Ne v]A14.3//m/[Ne ll]A12.8//m, and 2) [O lv]A25.9/im/[Ne ll]A12.8//m (Genzel et 
al. 1998), and b) optical: 1) [O m]A5007/H/3, 2) [N ll]A6583/Ha, 3) [S ll]AA6717,6731/H«, 
and 4) [O l]A6300/Ha (Sanders et al. 1988a; Armus et al. 1989; Veilleux et al. 1995; Due 
et al. 1997; Kim et al. 1998). As shown in Table 1, since the most MIR line ratios are 
upper-limit data, it is difficult to classify the ULIGs unambiguously solely with the ratios. 
Therefore, following Genzel et al. (1998), we use the diagram between the [O lv]/[Ne il] 
intensity ratio and the relative strength of 7.7 /im PAH feature (see Figure 5 in Genzel et 
al. 1998). This diagram shows that only two ULIGs, Mrk 231 and Mrk 273, are classified 
as AGN; i.e., higher than 50% of the energy comes from AGN. The optical classification 
is made with the criteria given in Figure 5 in Veilleux et al. (1995; see also Veilleux & 
Osterbrock 1987). 

The optical line ratios show that only four galaxies are classified as Hn-region (i.e., 
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starbursts) galaxies. Mrk 2310 and Mrk 273 are SI and S2, respectively. The remaining 
six ULIGs are LINERs. In Table 2, we compare the activity types based on the MIR 
classification with those based on the optical one. Both the MIR and optical classification 
give the same activity types for the following six galaxies; Mrk 231 and Mrk 273 (AGN), 
and IRAS 17208-0014, IRAS 20100-4156, IRAS 20551-4250, and IRAS 22491-1808 
(starburst) although IRAS 20551—4250 has been suspected to have LINER-like properties 
(Johansson 1991; Due et al. 1997). Since the MIR observations can probe the heavily- 
reddened, inner part of the nuclear region, we expect that the MIR spectroscopy would 
reveal hidden AGN in some ULIGs which are classified as starburst galaxies based on the 
optical observations. However, no such case is found. Contrary to the above expectation, 
six LINER-like ULIGs based on the optical observations turn to be starburst-dominated 
galaxies based on the MIR spectroscopy. This suggests strongly that these ULIGs which 
are optically classified as LINERs are not genuine AGN but shock-heated galaxies. If this 
is the case, we can conclude that both the MIR and optical classifications give the same 
activity types for all the ULIGs analyzed in this study. 

We mention the difference between the MIR and optical classification schemes. Since 
optical emission lines arise mostly from low-ionization ions, the optical classification uses 
the low-ionization emission-line ratios (Veilleux & Osterbrock 1987; Veilleux et al. 1995). 
On the other hand, Genzel et al. (1998) adopted the high-ionization lines such as [O iv] 
and [Ne v] in order to obtain firmer evidence for the central engine of AGN. Therefore, if 

3 The optical spectrum of Mrk 231 is dominated by very strong Fe II emission features, 
which are one of important characteristics of type 1 AGN (e.g., Boroson & Green 1992). 
Further, the presence of the broad Balmer emission lines provides the evidence for AGN 
in this galaxy (Boksenberg et al. 1977; Lipari, Colina, & Macchetto 1994 and references 
therein). 
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we would like to identify LINERs in future MIR spectroscopy, we will have to choose some 
low-ionization MIR emission lines such as [0 i]A63//m, [0 i]A145//m, and [C ll]A158/im 
(Spinoglio & Malkan 1992; Luhman et al. 1998). 

4. DISCUSSION 

The most important implication suggested by the MIR spectroscopy is that the 
majority of ULIGs have little evidence for AGN with high-ionization emission lines. This 
casts a serious problem; "Are low-ionization AGN classified with the optical spectroscopy 
really AGN ?" Since the discovery of LINERs (Heckman 1980), this question has been 
addressed many times (Heckman 1986; Ho 1998 and references therein). Here we consider 
what the case for the ULIGs isQ. There are three alternative ideas to explain the LINER-like 
spectra in the optical; a) a genuine AGN with low-ionization conditions (e.g., Ferland & 
Netzer 1983; Halpern & Steiner 1983), b) dense gas clouds photoionized by hot O stars 
(T eff £ 45,000 K) (Filippenko & Terlevich 1992; Shields 1992), and c) shock heating 
(Johansson & Bergvall 1985, 1988; Johansson 1991; Veilleux et al. 1995, 1997; Kim et al. 
1998; see also Heckman 1980; Dopita & Sutherland 1995). Since the MIR classification 
made by Genzel et al. (1998) deals only with high-ionization lines, it is impossible to probe 
genuine LINERs in the ULIGs. However, the most important point clarified by Genzel et al. 
(1998) is that the major energy source in the ULIGs is massive stars. This suggests strongly 
that the optical emission-line nebulae cannot be attributed mainly to the photoionization 
by the central engine of AGN. Therefore, we discuss the latter two possibilities. 

4 We do not discuss about LINERs associated with nuclei of ordinary galaxies in this 
paper. Those who are interested in this issue are recommended to see Ho, Filippenko, & 
Sargent (1997), Maoz et al. (1998), and Barth et al. (1998). 
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The intense starbursts are actually occurring in the heart of ULIGs (Armus et al. 
1989; Shaya et al. 1994; Scoville et al. 1998; Taniguchi et al. 1998). Further, the gas 
densities in them are generally higher than those in typical star forming regions in galaxies 
(Taniguchi & Shioya 1998). Therefore, we discuss the second possibility at first. Although 
LINERs are often regarded as low-luminosity extensions of typical luminous AGN, dense 
gas clouds photoionized by hot O stars (T c g ^ 45,000 K) also show LINER-like optical 
spectra (Filippenko & Terlevich 1992; Shields 1992). One of the important characteristics 
of such O-star LINERs is that the [O l]A6300 emission is significantly weak with respect 
to the Ea emission; [O I ]/Hct < 1/6. We mention that all the ULIGs except NGC 6240 
satisfy this condition (see Table 1). Although there is no firm evidence for hot O stars 
in the ULIGs, the presence of Wolf-Rayet stars is reported for some luminous infrared 
galaxies; e.g., IRAS 01003-2238 (Armus, Heckman, & Miley 1988) and NGC 1614 (Conti 
1991). If these high-temperature descendents of massive stars dominate the photoionization 
of ULIGs, they would lead to the LINER-like spectra. In Figure 1, we compare the optical 
emission-line ratios with the hot O-star photoionization models of both Filippenko & 
Terlevich (1992) and Shields (1992). The comparison shows that both the models cannot 
explain the data points of the LINER-like ULIGs. Although these models may explain the 
excitation properties of some starburst-dominated ULIGs, we conclude that hot O stars are 
not the major ionization source in the LINER-like ULIGs studied here. 

We consider the possibility of shock heating. It is naturally expected that some ULIGs 
are dominated by shock heating because they show evidence for the superwind activity; 
i.e., a blast wave driven by a collective effect of a large number of supernovae in the central 
region of galaxy mergers (Heckman, Armus, & Miley 1987, 1990; Veilleux et al. 1995; 
Heckman et al. 1996; Ohyama, Taniguchi, & Terlevich 1997; Kim et al. 1998). Further, very 
recently, Lutz et al. (1998) detected faint [O lv]A25.9/xm emission in a number of nearby 
starburst galaxies such as M82. This provides strong evidence for shock heating in the 
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starburst galaxies because they showed that the fast shock models by Dopita & Sutherland 
(1995) explain both the optical and MIR spectroscopic properties of M82 consistently. In 
order to examine if the shock heating models by Dopita & Sutherland (1995) can explain 
the optical emission-line properties of the ULIGs studied here, we compare the observed 
optical emission-line ratios with results of their shock models in Figure 1. It is shown 
that the majority of the ULIGs can be explained by the shock heating models with shock 
velocity between 100 km s _1 and 500 km s _1 (see also Taniguchi & Ohyama 1998). 

In summary, comparing the MIR and optical classifications for the sample of ULIGs, 
we have confirmed that the LINER-like ULIGs are not genuine AGN but shock-heated 
galaxies as suggested by Veilleux et al. (1995) and Kim et al. (1998). The recommendation 
should be to use both the optical and MIR classifications rather than one method. Finally, 
we should mention that our discussion is applicable only to infrared selected emission-line 
galaxies. Optically-selected LINERs may be powered by another energy source, e.g., AGN 
(see Ho 1998 and references therein). 
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Fig. 1— Excitation diagrams for the ULIGs except Mrk 231; a) [N Ilj/Ha (left), b) 
[S ll](A6716+6731)/Ha (middle), and c) [O l]A6300/Ha (right) are plotted against [O 
lll]A5007/H/3. Dashed curve shows the distinction among H II region-like objects, Seyfert 2 
(S2), and LINERs in each panel, taken from Veilleux et al. (1995). The alphabets correspond 
to those in Table 1. The hot O-star photoionization models by Filippenko & Terlevich (1992) 
are shown by dash-double dotted curves; for the solar abundances (thin) and for the 1.4 times 
solar abundances except nitrogen which is 4 times solar (thick). A logarithmic value of the 
ionization parameter is labeled at each asterisk mark. The composite hot O-star model by 
Shields (1992) is shown by a open square in each panel. The shock models are taken from 
Dopita & Sutherland (1995); 1) "shock + precursor" models (solid curves) with magnetic 
parameter B /n 1 / 2 = 0, 2, and 4 /xG cm 3 / 2 and shock velocity from 200 km s _1 to 500 km 
s _1 , and 2) "shock only" models (dash-dot curves) with magnetic parameter B/n 1 ^ 2 = 0, 2, 
and 4 /iG cm 3 / 2 and shock velocity from 150 km s _1 to 500 km s _1 . 
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Mrk 231 


< -0.22 


< -0.48 


A 










Sl s 


d 


Mrk 273 


-0.07 


-0.57 


A 


0.64 


-0.06 


-0.29 


-0.90 


S2 h 


e 


14348-1447 


< -0.68 


< -0.39 


H 


0.43 


-0.25 


-0.42 


-1.09 


L h 


f 


Arp 220 


< -1.51 


< -1.22 


H 


0.20 


0.36 


-0.13 


-0.79 


L f 


g 


NGC 6240 


-0.80 


< -1.30 


H 


0.07 


0.09 


0.04 


-0.43 


L c 


h 


17208-0014 


< -0.85 


< -0.72 


H 


-0.14 


-0.41 


-0.64 


-1.30 


H c,i 


i 


20100-4156 


< -0.21 


< -0.11 


H 


0.17 


-0.36 


-0.77 


-1.21 


H j 


j 


20551-4250 


< -0.92 


< -0.74 


H 


0.03 


-0.40 


-0.41 


-1.13 


H c,k 


k 


22491-1808 


< -0.35 


< -0.49 


H 


-0.19 


-0.37 


-0.62 


-1.24 


H c 


1 


23128-5919 


-1.43 


< -1.15 


H 


0.43 


-0.37 


-0.34 


-1.03 


L l,m 



a This ID is used to show each galaxy in Figure 1. 
b Logarithmic MIR emission-line ratios. 
c Logarithmic optical emission-line ratios. 

d A = AGN or AGN-like, SI = Seyfert 1, S2 = Seyfert 2, L = LINER, and H = H II or starburst galaxies. 
c Veilleux et al. (1995). 
f Armus et al. (1989). 

g Boksenberg et al. (1977), Sanders et al. (1988a), and Li'pari et al. (1994). 
h Sanders et al. (1988a). 

'Due et al. (1997) classified this galaxy as a LINER. 
JDuc et al. (1997). 

k Duc et al. (1997) and Johansson (1991) classified this galaxy as a LINER. 
'Young et al. (1996). See also Due et al. (1997). 

m Duc et al. (1997) classified this galaxy as a mixture of LINER/S2/starburst. 
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Table 2 

mir versus optical classification of the uligs 



MIR 


Optical 


ULIGs 


AGN 


Seyferts 


Mrk 231, Mrk 273 


AGN 


LINERs 




AGN 


H ii 




H ii 


Seyferts 




H ii 


LINERs 


Arp 220, NGC 6240, UGC 5101, 






12112+0305, 14348-1447, 23128-5919 


H ii 


H ii 


17208-0014, 20100-4156, 20551-4250, 22491-1808 
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